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1.0  EXECUTIVE  SUMMARY 


This  report  covers  work  perforaed  on  Phase  I  of  a  Small  Business 
Innovation  Research  project  for  AFOSRf  on  the  topic  of  cluster  beam  deposition 
of  thin  films.  The  project  was  intended  from  the  outset  to  provide  a  much 
firmer  scientific  basis  for  cluster  beam  deposition  than  has  existed  to  date. 
Weaknesses  which  were  apparent  in  the  literature  included:  1)  the  lack  of 
realistic  fluid  mechanical  treatments  in  models  of  the  cluster  growth 
environment;  and  2)  inadequate  methods  of  beam  characterlzat. ‘m.  It  was  felt 
that  Phase  I  work  could  result  in  significant  contributions  by  addressing 
these  problem  areas,  due  to  the  current  Interest  in  clusters  and  the  unique 
film  structures  which  may  be  possible  using  cluster  beams. 

Modeling  work  performed  on  the  project  used  the  method  of  characteristics 
to  generate  beam  expansion  rates,  including  some  off-axis  cases.  These  rates 
were  then  used  in  homogeneous  nucleatlon  and  growth  calculations.  It  was 
found  that  nucleatlon  rates  were  high,  but  that  clusters  could  grow  very 
little,  due  in  part  to  the  increase  in  temperature  experienced  through 
condensation  heating.  Beam  expansions  calculated  from  the  characteristics 
method  generally  led  to  somewhat  less  cluster  growth  than  when  simple  30** 
expansions  were  assumed,  as  has  been  common  in  the  literature. 

With  condensation  heating  neglected,  clusters  of  up  to  approximately  42 
atoms  could  be  predicted.  Preliminary  calculations  of  radiant  cooling  by  the 
clusters  were  performed  which  showed  that  a  potentially  Important  mechanism 
has  not  been  accounted  for  previously.  A  full  treatment  of  this  mechanism 
could  not  be  carried  out  within  the  constraints  of  the  project  but  may  be 
treated  in  a  subsequent  publication. 

A  general  conclusion  of  the  modeling  task  was  that  homogeneous  nucleatlon 
can  probably  condense  no  more  than  about  IX  of  the  metal  vapor  flux  into 
clusters.  This  conclusion  is  based  on  calculations  for  the  simple  nozzle  beam 
sources  described  in  the  literature,  and  for  the  time  being  neglects  the 
possibility  of  radiant  cooling. 

The  experimental  work  included  the  construction  and  testing  of  a  complete 
cluster  beam  apparatus,  including  nozzled  crucibles  and  heaters,  ionizing 
stage,  electrostatic  energy  analyzer,  and  ion  collector.  The  detection 
efficiency  of  the  apparatus  was  determined  for  metal  and  residual  gas  detection 
efficiencv.  By  such  methods  it  was  determined  that  a  beam  mass  fraction  of 
3.7  X  10"^  could  be  detected  for  all  clusters  in  the  size  range  of  400-600 
atoms  per  cluster,  and  5.6  x  10"^  for  all  clusters  in  the  800-1200  range. 

For  monoenergetic/monodlsperse  clusters,  the  limits  are  lower  by  about  a 
factor  of  ten. 

A  large  number  of  experiments  were  run  with  silver  and  aluminum  throughout 
the  vapor  pressure  regime  of  1  to  16  torr.  No  peaks  definitely  attributable 
to  clusters  were  seen  in  the  analyzer  data,  which  would  have  shown  clusters  of 
up  to  2000  atoms  per  cluster.  Thus,  conditions  for  cluster  production  and 
measurement  were  carefully  controlled;  if  upwards  to  30X  of  the  beam  can  be 
made  into  clusters  as  has  been  claimed  elsewhere,  they  should  have  been 
detectable. 
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In  spite  of  this  result »  It  was  found  that  the  noszled  source  did  produce 
a  considerable  aaount  of  snail  particles  which  could  be  seen  directly  or  with 
a  nlcroscope.  The  conditions  under  which  these  particles  were  produced  were 
studied.  These  studies  support  a  new  theory  of  heterogeneous  cluster 
fomatlon  on  Inner  crucible  surfaces  developed  by  Knauer.  The  heterogeneous 
theory  places  nore  stringent  requirements  on  crucible  temperature  control  than 
tile  homogeneous  mechanism;  It  Is  possible  that  particles  or  clusters  in  the 
desired  size  range  could  be  made  In  a  source  with  two  or  more  Independent 
temperature  zones. 

Workers  at  Hughes  have  recently  seen  evidence  of  600  atom  (3  nm)  clusters 
but  at  the  IZ  level  or  less.  One  option  for  future  cluster  deposition  work  Is 
to  attempt  a  more  rigorous  crucible  re*^eslgn,  In  order  to  optimize  for  the 
heterogeneous  cluster  growth  mechanism.  Two  factors  persuaded  us  not  to  take 
this  approach  further  than  the  effort  already  described  In  Section  5.4.2:  1) 

the  heterogeneous  theory  cannot  yet  provide  rigorous  design  Information,  and 
2)  Intuitively,  It  appears  unlikely  that  clusters  formed  on  surfaces  and  then 
entrained  In  the  nozzle  flow  could  comprise  a  significant  fraction  of  the  beam 
mass. 


An  alternate  approach,  based  on  growth  and  transport  of  ultraflne 
particles  or  clusters  In  a  quenching  gas,  appears  more  promising  In  our  view 
and  has  been  proposed  for  P^se  II  development.  This  report  described  a  test 
which  showed  that  such  particles  could  be  produced  In  our  apparatus.  The 
problems  of  gas/partlcle  transport  and  separation  appear  to  be  solvable  and 
are  addressed  In  the  Phase  II  proposal.  Most  of  the  theoretical  and 
experimental  capabilities  developed  In  Phase  I  work  can  be  applied  directly, 
since  nozzle  gas  flows  and  beam  characterization  will  continue  to  be  of 
central  Importance.  The  great  range  of  particle  sizes  obtainable  with 
gas**quench  techniques,  and  the  high  mass  utilization  factor,  should 
significantly  extend  the  possibilities  for  cluster  deposition  In  the  future. 


2.0  PHASE  I  WORK  STATEMENT 

The  Phase  I  Work  Statement  of  research  tasks  Is  shown  below  for  reference, 
as  given  In  the  contract  Instrument: 

(1)  Numerically  model  cluster  formation  using  nucleatlon  theory  or  the 
appropriate  kinetic  theory. 

(2)  Seek  experimental  confirmation  of  cluster  formation  In  collaboration. 

(3)  Ascertain  the  charged  species  present  In  the  Ionized  cluster  beam. 

(4)  Attempt  to  deposit  and  characterize  films  using  this  technique. 
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3.0  BACKGROUND  OF  IONIZED  CLUSTER  BEAM  DEPOSITION 


The  field  of  cluster  science  is  by  now  a  broad  one  and  no  atteapt  will  be 
■ade  to  review  it  here.  However*  it  is  useful  to  ezanlne  the  context  in  which 
cluster  bean  deposition  developed  and  to  conment  on  the  current  status  of  the 
technology. 


Cluster  science  began  as  a  fascination  with  things  that  happen  in 
supersonic  nossle  sources,  first  described  by  Kantrowitz  and  Gray  (1951). 

They  predicted  that  a  gas  passing  through  a  saall  orifice  into  vacuum  would 
cool  substantially,  transferring  aost  of  its  random  thermal  energy  into 
forward-directed  kinetic  energy.  This  led  to  atonic  and  molecular  beam 
sources  of  high  intensity.  It  was  later  found  that  under  certain  conditions  a 
highly  supersaturated  state  existed  in  the  beams,  leading  to  nucleation  and 
growth  of  gas  clusters.  Typical  sizes  ranged  from  10^  to  10^  molecules 
per  cluster.  Characterization  of  these  cluster  beams  became  an  active  area  of 
research  through  the  sixties  and  early  seventies.  Probably  the  best  review  of 
the  work  through  this  period  is  an  article  by  O.F.  Hagena  (1974). 


The  early  and  mid-seventies  also  saw  the  beginnings  of  a  strong  interest 
in  new  thin  film  technologies,  particularly  those  Involving  ion-activated 
deposition.  Taken  as  a  class,  these  methods  offered  better  film  strength, 
adhesion,  and  density.  Certain  drawbacks  were  also  noted  such  as  radiation 
damage  to  electronic  circuits  during  metallization  processes,  implantation  of 
ionized  impurities  into  films,  and  difficulties  in  controlling  beams  with  low 
mass-to-charge  ratios. 


Work  on  Ionized  Cluster  Beam  deposition  was  begun  in  this  period  by  Takagl 
and  co-workers  at  Kyoto  University  (Takagl,  1975).  A  nozzle  source  similar  to 
those  used  for  gas  clusters  was  found  to  produce  large  metal  clusters  as 
well.  By  charging  massive  clusters  instead  of  individual  atoms,  it  was 
reasoned,  one  could  obtain  activated  deposition  without  entailing  any  of  the 
drawbacks  noted  above.  A  typical  source  apparatus  including  cluster 
ionization  and  acceleration  components  is  shown  in  Figure  1.  A  number  of 
application-oriented  papers  by  Takagl 's  group  showed  considerable  promise  for 
the  technique,  even  though  the  distinguishing  mechanism  for  film  growth  had 
not  been  conclusively  demonstrated.  The  mechanism  proposed  by  Takagl  was  that 
of  total  breakup  of  the  clusters  on  impact  and  partitioning  of  the  cluster 
energy  into  surface  migration  energy  of  atoms,  as  shown  in  Figure  2. 

Another  mode  of  growth  involving  compaction  and  flattening  has  recently 
been  predicted  through  molecular  dynamics  calculations  (Mullev ,  1987) . 

Examples  of  his  results  are  shown  in  Figure  3.  The  molecular  dynamics 
approach  measures  surface  mobility  of  individual  atoms,  but  it  is  interesting 
to  observe  that  something  like  a  plasma  sprayed  coating  might  be  possible  on  a 
nanometer  scale.  This  could  lead  to  better  films  by  reducing  the  tendency  to 
form  loose  columns  during  growth. 

The  results  reported  by  Takagl  began  to  arouse  Interest  outside  Japan  in 
the  early  eighties.  ICB  systems  were  being  used  on  a  commercial  scale  in 
Japan  and  a  licensing  agreement  with  the  Eaton  Corporation  made  a  system 
available  in  the  U.S.  Unfortunately,  the  Eaton  systems  lacked  the  sort  of 
beam  characterization  facilities  which  most  researchers  required.  It  became 
clear  that  laboratories  interested  in  the  potential  benefits  of  cluster 
deposition  would  have  to  develop  their  own  techniques  and  equipment. 
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Figure  1.  T3rplcal  apparatus  for  ionized 
cluster  beam  deposition. 
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Figure  2.  Film  growth  by  cluster  fragmentation. 


Figure  3.  Film  growth  by  cluster  deformation  (a)  with 
E"0.2  eV/atom,  (b)  with  E«3  eV/atom  (From  Huller,  1987). 
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since  then  there  has  also  been  some  controversy  regarding  the  possibility 
of  condensing  metals  In  a  pure  vapor  expansion.  The  primary  difficulty  Is 
that  a  great  deal  of  latent  heat  must  be  released  In  forming  clusters.  In 
most  cluster  sources  this  heat  Is  gradually  taken  up  by  a  large  number  of 
Inelastic  collisions.  A  pure  vapor  can  take  up  some  of  Its  own  latent  heat  In 
this  way»  but  It  Is  far  easier  to  "seed”  metal  clusters  In  a  flow  of  Inert 
quenching  gas.  Sources  which  use  an  oven  or  pulsed  laser  to  vaporize  a  metal 
Into  a  gas  are  now  a  mainstay  of  research  In  the  chemistry  and  physics  of 
small  clusters. 

Takagl  has  steered  away  from  the  gas-quenchi'^g  approach,  possibly  to  avoid 
the  apparent  complication  due  to  gas  flows.  A  number  of  phenomena  such  as 
convection  and  nozzle  clogging  by  cluster  Impaction  must  be  considered. 
However,  these  are  essentially  transport  problems.  The  rellabllty  and 
efficiency  with  which  gas-quenched  sources  produce  clusters  leads  us  to 
believe  that  ways  should  be  found  to  exploit  these  as  well.  This  has  been 
tried  with  some  success  by  Kulper  (1981).  We  have  proposed  for  Phase  II  work 
a  variant  on  this  technique,  employing  a  more  sophisticated  approach  to  the 
problems  of  gas  flow  and  control.  Because  larger  clusters  or  particles 
(greater  than  5  nm)  can  also  be  made  by  such  methods,  the  deformation  mode  of 
film  growth  should  apply  (Figure  2),  thus  offering  more  radical  departures 
from  current  deposition  technologies. 


A.O  MODELING  OF  NOZZLE  BEAM  EXPANSION  AND  CLUSTER  GROWTH 
4.1.  Nucleatlon  and  Growth  Model 

We  have  performed  numerical  modeling  of  the  physical  phenomena  involved  in 
the  process  of  nucleatlon  when  a  gas  expands  Into  vacuum  through  a  supersonic 
nozzle.  For  a  gas  to  condense,  it  has  to  be  brought  to  a  supersaturation 
state  where  the  pressure  and  temperature  of  the  gas  fall  into  the  region  of 
the  P-T  plane  where  only  liquid  phase  may  exist  under  equilibrium  conditions. 
This  generally  leads  to  some  combination  of  gas  cooling  and  a  rise  in  pressure 
to  bring  the  system  toward  an  equilibrium  state.  In  the  case  of  a  one 
dimensional  gas  flow  through  a  supersonic  nozzle,  the  thermodynamic  state  of 
the  gas  Is  determined  entirely  by  the  cross-sectional  area  of  the  beam  If  the 
expansion  may  be  considered  Isotropic.  The  cross-sectional  area  A(x)  Is 
related  to  the  Mach  number  M(x)  by  the  equation 

_ iJi- 

1  r  2  vl  2"! 

A(x)/A(0)  -  Mt^)  I  ^  (1  +  ^  )J  (1) 

where  x  Is  the  axial  distance  from  the  nozzle,  A(0)  Is  the  cross-sectional 
area  of  the  nozzle  throat,  and  y  is  the  specific  heat  ratio.  From  equation 
(1),  the  Mach  number  at  the  throat  Is  M(0)-1.  M(x)  can  be  computed  If  the 
area  ratio  A(x)/A(0)  Is  known  as  a  function  of  the  distance  x  from  the 
throat.  The  pressure  and  the  temperature  of  the  gas  are  then  determined  with 
the  following  gasdynamlcs  formulas: 
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r  5  -|-y/(y"1) 

p  -  1 1  +  (•ri)M^/2j 

(2) 

T  -  [l  +  (r-DM^/a] 

(3) 

where  P^  and  are  stagnation  pressure  and  temperature,  l.e.,  pressure 

and  temperature  at  the  source  where  the  gas  velocity  Is  zero.  For  the  flow 
beyond  the  sonic  point  (x“0)  the  Mach  number  is  typically  greater  than  unity 
and  increases  rapidly  with  the  distance  from  the  throat,  equation  (3).  The 
gas  Is  then  cooled  rapidly  and  may  reach  a  supersaturated  state.  Such 
qualitative  reasoning  forms  the  basis  of  the  explanation  of  the  experimentally 
obsetrved  phenomenon  of  cluster  formation  In  nozzle  beam  sources. 

To  describe  the  process  of  cluster  formation  In  a  more  quantitative  way, 
we  have  adopted  the  following  approach.  The  expansion  of  the  gas  Is  assumed 
to  be  a  one-dlmenslonal  steady  state  flow  so  that  all  physical  quantities 
depend  only  on  the  variable  x  which  measures  distance  from  the  nozzle  throat . 
The  governing  equations  for  the  one-dlmenslonal  flow  Include  the  continuity, 
momentum,  and  energy  equations  and  contain  the  effects  of  nucleatlon  as 
obtained  from  the  classical  nucleatlon  theory,  as  well  as  cluster  growth  due 
to  condensation.  These  equations  are  (Hill,  1966): 

^ 

Y  YM  J 

(4) 

".t[x 

(5) 

(6) 

(7) 

gl-  'o' 

(8) 

“  ^3  ^  ^  ^  'o 

(9) 

||3  -  Sir  I  A  (10) 
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In  the  above  equations,  P  is  pressure,  T  the  teaperature,  M  the  Mach 
nuaber,  y  is  the  ratio  by  aass  of  liquid  and  vapor  phases,  A  the  beaa 
cross-sectional  area,  X*h/Cp  T  where  h  is  the  enthalpy  of  the  gas-liquid 
transition  and  Cp  is  the  specific  heat  at  constant  pressure,  y  is  the 
ratio  of  specific  heats  at  constant  pressure  and  voluae,  is  the  density 
of  the  liquid,  a  is  the  total  flow  rate,  I  is  the  nucleation  rate,  r^  is  the 
radius  of  the  nucleating  drops,  u  is  the  velocity  of  the  gas,  and  dr/dt  is  the 
drop  radius  growth  rate.  Q]^,  Q2>  and  Q3  are  auxiliary  variables  related 
to  the  total  nuaber  of  condensed  droplets  and  their  total  surface  area  so  that 
the  surface  averaged  drop  radius  is 


T  - 


(11) 


We  assuaed  that  at  location  z,  droplets  with  the  radius  rQ(z)  are 
spontaneously  created  at  the  rate  (per  unit  voluae)  I(z)  given  by  the 
classical  nucleation  theory  (Frenkel,  1946): 

I  ■  /  2af  m  exp  ^-Airo^r  /3kT^  (12) 

where  a  is  the  surface  tension  of  the  liquid,  a  is  the  aass  of  one  aolecule, 
V  is  the  voluae  of  one  aolecule  in  the  liquid,  k  is  the  Boltzaann  constant, 
and  Tf>  is  the  critical  radius  at  teaperature  T  and  pressure  P.  The  critical 
radius  is  given  by 


'c  "  f^^&TSrPTTg) 

where  Pg(T)  is  the  flat-film  saturation  pressure.  The  drop  growth  rate  is 
then  determined  from  the  equations  relating  the  net  mass  flow  to  and  from  the 
droplet  surface  due  to  condensation  and  re-evaporation. 


dr 

dt 


(14) 


where  ^  is  the  condensation  coefficient  (usually  assumed  to  be  close  to 
unity),  and  6  and  6p  are  the  condensation  and  evaporation  rates, 
respectively.  The  condensation  rate  is  the  mass  flux  to  the  surface  of  the 
drop  (per  unit  area,  per  unit  time)  and  is  given  by  the  expression  derived 
from  the  kinetic  theory 


6  -  P/  ✓'(2it  k  T/m) 

e 

Similarly,  the  evaporation  rate  is  given  by 


(15) 


Bp  -  Pp/  /  (2Tr  k  Tp/m) 


(16) 
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vbere  Pq  and  Tj)  refer  to  the  pressure  of  the  liquid  drop  and  to  Its 
temperature y  respectively.  The  pressure  Pq  Is  the  pressure  of  the 
equilibrium  vapor-liquid  state  at  temperature  Tp  corrected  for  the  effects 
of  the  surface  tension  for  the  drop  with  radius  r: 

Pd  "  Ps(T>  exp  (2am/(pLkTDr)).  (17) 


The  temperature  of  the  drop  Is  usually  higher  than  the  temperature  of  the 
condensing  vapor  because  of  the  latent  heat  release.  7%is  can  amount  to 
differences  of  tens  or  hundreds  of  degrees,  depending  on  the  latent  heat  of 
condensation,  and  thus  strongly  affects  the  re-evaporatlon  rate.  This  effect 
has  not  been  accounted  for  In  some  of  the  recent  literature  on  cluster  beam 
modeling.  This  points  out  the  vlrttial  necessity  of  removing  the  latent  heat 
by  collisions  with  an  Inert  background  gas,  as  we  plan  to  do  in  the  Phase  II 
work.  The  drop  temperature  Is  computed  from  the  following  algebraic  equation: 


KC(l-ajjTjj/ST)  -  K(l-0  a  -l)+  C(l-Sp/e)  (X-1)^  -0  (18) 

where  K  ■  ( yH)/2( y-l)  is  the  average  specific  energy  of  the  arriving 
molecules  and  oi  Is  the  thermal  accommodation  coefficient  which  Is  always 
close  to  the  unity. 


Finally,  it  Is  assumed  that  the  vapor  satisfies  the  perfect  gas  law 
P  ■  p  kT/m 


(19) 


which  also  gives  a  simple  expression  for  the  speed  of  sound 
c  ■  /  y  kT/m 

and  the  Mach  number  M“u/c. 


(20) 


The  Numerical  Method 

Equations  (4)-(10)  with  (16)  and  (18)  describe  one  dimensional  gas  flow 
with  condensation.  To  solve  these  equations,  numerical  methods  must  be  used. 
Equation  (18)  was  solved  by  the  bisection  and  secant  method  (Oahlqulst  & 

Bjork,  1974)  and  equations  (4)-(l0)  by  the  fourth  order  Runge-Eutta  method,  as 
•  suggested  by  Hill  (1966). 

The  numerical  method  was  tested  by  repeating  some  of  the  calculations 
reported  by  Hill  (1966)  for  water  vapor  condensing  In  a  supersonic  nozzle  of 
prescribed  geometry.  Agreement  was  good,  predicting  droplet  formation  of  3.7X 
of  the  vapor  with  average  droplet  size  of  17  nm.  These  numbers  correlated 
A  well  with  experiments  referenced  by  Hill.  Area  distributions  for  different 
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nocsles  considered  by  Hill  are  reproduced  here  as  Figure  4.  Figure  5  shows 
results  of  his  calculations  for  different  values  of  the  stagnation  pressure 
and  temperature  for  one  of  those  nossles.  In  Figures  6  and  7,  we  show  results 
of  our  calculations  for  two  cases  denoted  as  1  and  4  In  Figure  5.  Both 
figures  clearly  show  departures  of  the  computed  pressures  from  the  Isentroplc 
(lower)  curves.  Indicating  the  onset  of  condensation.  Our  calculations  are 
also  in  good  agreement  with  the  results  of  Rill  (Figure  5).  The  only 
difference  Is  that  In  our  calculations,  the  condensation  seems  to  start 
slightly  earlier.  The  probable  reason  for  this  discrepancy  Is  our 
approximation  of  the  area  distribution  of  the  experimental  nozzle  from 
Figure  4  by  a  parabolic  expansion  fitted  at  x~0.5  and  10cm,  which  leads  to  a 
slightly  higher  expansion  ratio  than  the  original  nozzle.  Considering  other 
sources  of  discrepancy  such  as  the  possible  use  of  different  expressions  for 
the  saturation  pressure  (not  given  by  Hill)  In  the  literature,  and  scatter  In 
the  experimental  data,  the  agreement  between  the  two  sets  of  calculations  Is 
sufficient  to  validate  our  numerical  code. 

4.2  Vapor  Flow  Model 


The  next  step  was  to  properly  calculate  the  rate  of  expansion  of  a  free 
vapor  Jet,  an  approach  not  yet  taken  elsewhere.  The  expansion  of  a  molecular 
beam  Into  vacuum  Is  very  different  from  the  flow  through  a  long  diverging 
nozzle  since  there  are  no  walls  to  keep  the  flow  within  prescribed 
boundaries.  This  difficulty  Is  usually  circumvented  by  assuming  that  the 
expanding  flow  has  a  conical  shape  %rlthln  a  prescribed  apex  angle  of  about 
20-30°  (Usui  et  al.,  1985,  Yang  &  Lu,  1986).  In  this  research,  we  used  two 
different  ways  of  determining  the  area  distribution  of  the  expanding 
supersonic  flows.  One  Is  consistent  with  the  above  described  assumption  of 
the  conical  flow  and  the  other  determines  the  area  distribution  from  equation 
(1)  with  the  Mach  number  obtained  from  the  solution  of  the  full  gasdynamlcs 
equations  for  the  axlsymmetrlc ,  Isentroplc  flow  expanding  from  the  cylindrical 
nozzle  Into  vacuum.  The  latter  approach,  which  should  be  more  realistic,  will 
be  briefly  described  below. 

The  flow  of  the  metallic  vapors  through  the  orifice  was  modeled  as  a 
two-dimensional  axlsymmetrlc  nozzle  flow  using  the  method  of  characteristics. 
This  method  provides  detailed  flow  Information,  Including  the  angular 
distribution  of  mass  flux.  The  calculation  begins  at  the  most  constricted 
part  of  the  flow  (the  nozzle  exit)  which  Is  the  sonic  point  (M>1).  Because 
the  gas  flow  beyond  this  point  Is  decoupled  from  disturbances  upstream,  It  Is 
only  necessary  to  find  P  and  T  at  M*1  to  carry  out  the  computations  In  the 
expansion  region.  These  quantities  are  obtained  from  equations  (2)  and  (3) 
above.  While  Is  Is  possible  that  some  nucleatlon  and  growth  occurs  before  the 
sonic  point,  calculations  by  Yang  and  Lu  (1986)  Indicated  that  subsonic 
clustering  was  almost  negligible. 
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p/p. 


Figure  4.  Area  expansion  ratios  for  diverging 
nozzles  in  water  vapor  experlaents  (fron  Hill,  1966). 


Figure  5.  Experlaental  data  of  Binnie  and  Green  eoapared 
with  nucleation  and  growth  calculations  aade  by  Hill  (1966) 
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Figure  7.  Valldatloa  of  cluater  growth  code,  ualng  conditions 
for  Case  4  of  Figure  5  (Blnnle  and  Green  notsle, 
pQ  -  0.646  at».  To  -  127*C). 
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The  method  of  characteristics  described  by  Shapiro  (1954)  la  used  for 
exact  stepwise  numerical  calculation  of  azlsymmetrlc  supersonic  flows.  The 
flow  Is  also  assumed  to  be  steady  and  Irrotatlonal.  The  method  of 
characteristics  solves  for  the  flow  field  parameters  by  following 
rlgiht'Tunnlng  and  left-running  Mach  lines  (characteristics)  which  trace 
azlsymmetrlc  compression  and  expansion  waves.  Because  the  waves  weaken  as 
they  move  away  from  the  axis  of  symmetry  and  strengthen  as  they  move  towards 
it*  the  equations  of  motion  must  be  Integrated  for  the  strength  of  the  waves 
as  well  as  their  location.  Thus  there  are  four  unknowns  to  solve  for  at  each 
intersection  of  right-running  and  left-running  characteristics:  the  axial 
location  X,  the  radial  location  r,  the  velocity  V  and  the  flow  directional 
angle  9.  The  four  equations  for  the  unknowns  are  written  In  terms  of  known 
Information  at  another  point  on  the  right-running  characteristic  and  another 
point  on  the  left-running  characteristic.  When  parameters  vary  along  the 
characteristic  between  the  known  point  and  the  unknown  point ,  the  average 
value  Is  used. 

Depending  upon  what  additional  Information  is  known  at  the  unknown  point , 
the  solution  of  the  four  equations  takes  various  forms: 

(1)  The  basic  case  Is  the  free  Intersection  of  two  characteristics. 

(2)  If  the  unknown  point  Is  on  the  axis  of  symmetry »  then  r  and  6  are  known 
to  be  zero,  so  the  problem  Is  reduced  to  two  unknowns. 

(3)  If  one  of  the  known  points  Is  on  the  axis  of  symmetry,  then  a  special  form 
of  the  equations  for  the  four  unknowns  must  be  used  because  the  equations  are 
Indeterminate  on  the  axis.  The  limit  of  the  Indeterminate  value  as  the 
characteristic  approaches  the  axis  is  used. 

(4)  If  the  unknown  point  Is  on  the  wall,  then  the  flow  angle  Is  known  from  the 
wall  angle.  Three  equations  are  solved  for  x,  r,  and  V. 

(5)  If  the  characteristic  runs  uninterrupted  from  the  axis  of  symmetry  to  the 
wall,  then  yet  another  form  of  the  equations  Is  used,  incorporating  the  limits 
on  the  axis  and  the  wall  angle. 

(6)  If  the  unknown  point  la  on  the  edge  of  the  flow  field,  where  In  actuality 
a  mixing  layer  exists,  then  the  flow  velocity  can  be  found  from  matching  the 
jet’s  pressure  with  the  ambient  pressure.  This  leaves  three  equations  for  x, 
r,  and  9. 

The  characteristics  method  is  useful  only  in  regions  where  the  mean  free 
path  between  collisions  Is  small.  The  mean  free  path  is  Inversely 
proportional  to  the  number  density  of  the  molecules  and  the  effective 
collision  cross  section.  In  the  vapor  above  liquid  aluminum  at  2139  R  and 
20  torr,  the  mean  free  path  Is  about  0.0012  cm.  When  the  vapor  Is  expanded 
through  a  nozzle,  however,  the  mean  free  path  increases  to  about  0.25  cm  at 
the  point  where  M“10.  When  the  mean  free  path  approaches  the  physical 
dimensions  of  the  device,  further  fluid  flow  calculations  become  meaningless, 
since  the  trajectories  of  the  molecules  will  be  unchanged  by  Intermolecular 
collisions  (a  condition  known  as  free  molecular  flow).  Since  the  method  of 
characteristics  experienced  convergence  stability  problems  at  high  Mach 
numbers,  the  calculation  was  terminated  when  the  Mach  number  exceeded  the 
ll*lt  for  free  molecular  flow.  This  generally  did  not  occur  until  the  flow 
was  many  nozzle  diameters  beyond  the  orifice. 
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Primarily •  the  method  of  characteristics  code  was  used  to  generate  the 
flow  expansion  ratio  A(x)/A(0)  along  streantubes  for  use  In  the  condensation 
computer  model  (as  Input  to  equation  (D).  The  flow  area  A(z)  of  the 
streamtube  Is  referenced  to  the  flow  area  of  the  streamtube  at  the  sonic  point 
A(0)»  where  the  flow  velocity  Is  equal  to  the  speed  of  sound.  Figure  8  shows 
several  streamtube  area  histories  from  the  2-0  method  of  characteristics 
calculation  compared  with  the  usual  choice  of  one-dlmenslonal  flow  In  a  30 
degree  cone.  The  2-D  area  histories  are  seen  to  depend  strongly  on  the 
starting  position  (radial  distance  from  centerline)  of  the  streamtube  In  the 
orifice.  This  Information  was  used  In  the  cluster  growth  code  to  determine 
whether  significant  differences  In  the  sizes  or  numbers  of  clusters  should 
occur  off  axis. 

The  final  flow  angle  of  the  streamtubes  Is  another  useful  output  of  the 
characteristics  model.  This  Information  can  be  used  to  estimate  the  mass 
deposition  rate  as  a  function  of  position  on  a  substrate.  The  thickness 
profile  of  a  deposited  film  can  then  be  used  to  test  the  model;  experimental 
data  showed  good  agreement  with  predictions,  as  will  be  shown  In  a  later 
section. 

The  mass  flux  of  the  whole  orifice  flow  field  Is  shown  In  Figure  9.  The 
length  of  the  vectors  Is  proportional  to  the  mass  flux  ■  pft  where  p  Is 
the  density  and  v  Is  velocity.  The  direction  of  the  vectors  Is  the  flow 
angle.  It  can  be  seen  that  the  orifice  flow  at  20  torr  is  peaked  toward  the 
center  but  is  not  strongly  directional. 

Expansion  ratios  were  calculated  for  ■  1,  2,  5,  10  and  20  torr,  but 
were  found  to  be  virtually  identical  for  axial  streamtubes.  In  reality,  some 
differences  should  appear  idien  Pq  Is  a  few  torr  or  less  because  of  boundary 
layer  effects  In  the  nozzle.  The  nucleatlon  and  growth  work  gave  significant 
condensation  only  at  the  higher  pressure,  however,  where  such  effects  should 
be  negligible. 

In  Figure  10,  we  present  results  for  the  area  expansion  ratio  for  aluminum 
vapor  expanding  Into  vacuum  through  a  cylindrical  nozzle  1  mm  In  diameter.  In 
this  figure,  results  computed  from  the  axisymmetrlc  method  of  characteristics 
are  compared  with  the  area  distribution  for  the  conical  flow  of  30®  apex 
angle.  It  can  be  seen  that  within  about  0.5mm  from  the  nozzle,  the  expansion 
rate  dA/dx  is  greater  for  the  conical  flow  than  for  the  free  flow.  Beyond 
that  point,  however,  the  expansion  ratio  of  the  free  flow  Is  substantially 
greater  than  for  the  conical  flow  for  all  stagnation  pressures  considered. 
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Figure  8.  Areal  expanaioo  of  etreaatubea  tb.  initial 
radial  atarting  diatance  froa  nocale  centerline. 
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Figure  9.  Distribution  of  ■ass  flux  in  free 
nozzle  expansion  froa  20  torr  crucible. 
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4.3  Results  of  Calculations 


We  chose  aluminum  as  a  working  material  In  our  numerical  simulations 
because  the  experimental  work  was  Initially  focused  on  aluminum  vapor  beams. 
Table  I  contains  our  primary  results.  In  this  Table  under  the  heading  "Flow" 
we  have:  flow  through  diverging  nozzles,  as  studied  by  Hill;  conical  flow  Into 
30^;  and  free  expanding  flow  calculated  by  the  method  of  characteristics. 

Pq  and  Tq  are  the  crucible  stagnation  pressure  and  temperature,  d  Is  the 
nozzle  diameter,  and  y  Is  the  mass  fraction  of  vapor  which  condenses  Into 
clusters.  The  last  two  columns,  rj[  and  rf ,  give  the  Initial  radius  of 
clusters  which  nucleate  at  the  throat  and  their  final  radius  after  reaching 
the  end  of  the  computational  zone,  respectively.  Besides  the  cases  listed, 
calculations  were  also  performed  for  aluminum  at  Pq  ■*  1,  2,  5,  and  10  torr, 
and  for  off-axis  streamtubes.  All  of  the  cases  shown  are  for  flow  near  the 
beam  axis. 

In  Figure  11,  we  present  results  for  temperature  (normalized  to  the 
stagnation  value)  as  a  function  of  the  distance  from  the  nozzle  for  the 
stagnation  pressure  Pq  ■  20  torr,  the  nozzle  diameter  d  ■  0.1  cm,  and  a 
conical  flow  of  30*^  apex  angle  (run  3  In  Table  I).  Here  the  lower  curve 
corresponds  to  the  Isentroplc  expansion  without  condensation  and  the  upper 
curve  Is  from  solution  of  equations  (4)-(10)  which  Include  the  effects  of 
condensation.  The  conical  flow  case  was  run  for  comparison  with  the  results 
of  Usui  et  al.  (1985)  who  assumed  a  simple  30°  expansion. 

The  fact  that  the  curves  In  Figure  11  do  not  coincide  Indicates  that  there 
Is  some  condensation  with  latent  heat  release  occurring.  This  fact  alone, 
however,  Is  not  sufficient  to  draw  the  conclusion  that  clusters  are  growing. 

In  Figure  12  we  show  the  radius  of  the  nucleating  clusters  at  position  x.  It 
Is  assumed  that  the  clusters  begin  with  radius  rg  (x)  •  1.3rc(x).  (The 
multiplier  of  1.3  corresponds  to  a  cluster  size  for  which  the  growth 
probability  Is  near  lOOZ.)  It  can  be  seen  that  the  critical  radius  decreases 
with  X  due  to  the  lower  vapor  temperature.  Also  shown  Is  a  curve  that  gives 
the  final  radius  of  the  clusters  that  originated  at  x  and  traveled  to  the  end 
of  the  computational  domain;  the  vertical  displacement  of  the  curves  Indicates 
the  growth  shown  by  a  given  cluster.  Very  clearly  these  results  show  that 
there  was  little  cluster  growth  In  run  3.  The  departure  of  the  curve  with 
condensation  In  Figure  11  from  the  isentroplc  one  Is  therefore  caused  by  the 
latent  heat  release  due  to  the  formation  of  a  very  large  number  of  small 
nuclei,  for  which  the  condensation  rate  Is  only  slightly  greater  than  the 
evaporation  rate.  A  similar  conclusion  is  reached  for  run  4,  where  a  larger 
orifice  was  assumed  (0.2  cm).  The  beam  temperature  profile  Is  shown  In  Figure 
13  and  droplet  growth  In  Figure  14.  The  mass  fraction  y  Is  somewhat  higher 
as  shown  In  Table  I,  In  agreement  with  the  general  finding  that  high  values  of 
Pnd  lead  to  more  condensation  (Hagena,  1974).  We  should  note  that  Usui  et 
al.  (1985)  performed  similar  calculations  for  the  same  conditions,  but  In 
their  paper  the  curve  with  condensation  Is  removed  further  from  the  Isentroplc 
curve  than  In  our  case.  Their  results  were  thus  more  optimistic  than  ours. 

The  reasons  for  this  discrepancy  are  not  known  at  the  present  time,  since 
their  paper  states  only  final  results  without  any  discussion  of  the  basic 
equations  and  the  numerical  method  so  that  the  comparison  Is  difficult. 
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Figure  11.  Calculated  beaa  temperature  va.  distance 
froB  0.1  ca  orifice  nozzle.  For  aluainua  in  30^  expansion 
with  Pq  ■  20  torr,  Tg  ■  2140K.  Syabola  are  for  curve 
labeling  only:  squares  for  isentrope,  circles  for  condensing  vapor. 
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Figure  12.  Cluster  starting  radius  vs.  position  at  which  it 
nucleates  (lower  curve)  and  corresponding  final  radius 
(upper  curve).  For  30**  expansion  froa  ■  20  torr, 

0.1  cw  nozzle  diaaeter. 
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Figure  13.  Calculated  bean  tenperature  va.  diatance 
from  0.2  cm  orifice  nossle.  For  aluainua  in  30**  ezpanaion 
with  Pq  -  20  torr,  Tq  •  2140R.  Syabola  are  for  curve 
labeling  only:  squarea  for  iaentrope,  circles  for  condensing  vapor. 
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Figure  14.  Cluster  starting  rsdl\w  vs.  position  at  which  It 
nucleates  (lower  curve)  and  corresponding  final  radius 
(upper  curve).  For  30®  expansion  with  Pq  •  20  torr, 

0.2  cw  noscle  dlaweter. 
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Run  5  used  area  expansion  ratios*  calculated  from  the  characteristics 
model  rather  than  the  30°  assumption.  Temperature  and  cluster  size  are 
shown  In  Figures  15  and  16.  Again*  little  or  no  growth  occurred  even  though 
the  nucleatlon  rate  was  high.  However*  the  condensed  mass  fraction  dropped  by 
30Z  compared  with  case  3  (30°  expansion).  In  these  runs,  despite 
condensation  reaching  between  1  and  2Z  of  the  mass  of  the  vapor  (see  Table  1) 
the  radius  growth  Is  negligible.  These  cases  should  be  contrasted  with  runs  1 
and  2  for  water  vapor  where  the  radius  of  the  nucleating  droplets  grows  by  an 
order  of  magnitude  In  the  condensation  zone. 

The  relatively  small  radius  growth  rate  can  be  explained  on  the  basis  of 
the  following  reasoning.  An  upper  limit  to  the  growth  rate  Is  found  by 
setting  Rp  to  zero*  thus  neglecting  re^evaporatlon.  The  condensation  rate 
(15)  will  have  Its  maximum  value  just  beyond  the  throat  since  the  pressure 
decreases  more  rapidly  than  the  square  root  of  temperature.  For  Pq  ■  20 
torr  and  Tq  ••  2140°K  the  maximum  value  of  the  growth  rate  computed  for 
aluminum  Is  dr/dt*0.04  cm/s.  The  velocity  of  the  clusters  Is  comparable  to 
the  velocity  of  sound  at  the  throat  which  from  (20)  Is  c>9xl0^  cm/s. 

Therefore,  under  the  most  favorable  conditions*  one  may  expect  clusters  to 
grow  from  rj^  ■  0.3  nm  to  rf  ■  1.6  nm  (1000  A1  atoms)  In  a  distance  of  0.33 
cm.  In  reality,  dr/dt  will  decrease  with  distance  so  that  the  condensation 
zone  Is  always  much  less  than  1  cm;  the  maximum  growth  rate  cannot  be 
sustained  over  a  distance  much  greater  than  about  one  nozzle  diameter  (0.1  cm 
In  this  example).  Nucleating  droplets  probably  do  not  grow  by  more  than  about 
0.4  nm  even  with  evaporation  negelected*  yielding  clusters  of  perhaps  90  atoms. 

To  confirm  these  estimates,  we  have  rerun  cases  3-5  with  Bp  set  to  zero 
In  our  program  (runs  6-8).  Results  are  shown  in  Figures  16-22.  The  total 
amount  of  condensed  phase  Increases  If  evaporation  Is  neglected,  but  even  then 
the  radius  increase  does  not  exceed  0.3  nm.  To  be  more  exact,  the  largest 
radius  achieved  In  our  simulations  is  0.55  nm  which  corresponds  to  a  cluster 
consisting  of  42  aluminum  atoms  at  bulk  density.  This  Is  consistent  with 
experimental  results  reported  by  Kulper  (1981)  for  silver  beams. 

4.4  Limitations  of  Model 


Since  In  Takagl's  experiments  clusters  consisting  of  500-1000  molecules 
are  reported,  it  Is  worth  discussing  possible  limitations  In  our  modeling  that 
could  lead  to  overly  pessimistic  results.  Probably  the  weakest  point  In  our 
calculations  Is  the  assumption  that  the  vapor  nucleates  at  a  point  x  with  a 
monochromatic  spectrum  of  droplets  of  radius  1.3r(..  If  these  droplets  grow 
very  rapidly  as  In  the  case  of  water  vapor  In  a  supersonic  nozzle,  this 
assumption  Is  not  important  since  the  size  of  the  droplets  after  passing  the 
condensation  zone  Is  to  a  large  extent  independent  of  the  Initial  size  (before 
entering  the  condensation  zone).  In  reality,  however,  the  nucleatlon  process 
creates  an  entire  spectrum  of  drop  sizes  which  may  be  determined  from  the 
Gibbs  distribution.  If  they  grow  only  very  little,  the  final  distribution  of 
sizes  may  depend  strongly  on  the  Initial  distribution.  Moreover,  since  larger 
droplets  grow  faster  (the  growth  In  volume  Is  proportional  to  4mr^dr/dt) 
the  presence  of  larger  droplets  In  the  initial  distribution  could  lead  to  a 
prediction  of  much  larger  clusters. 
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Figure  15.  Calculated  bean  temperature  va.  diatance 
from  0.1  cm  orifice  nosale.  For  aluminum  by  method  of  characterietlcs 
with  Pq  ■  20  torr,  Tg  ■  2140R.  Symbols  are  for  curre 
labeling  only:  squares  for  Isentrope,  circles  for  condensing  vapor. 
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Figure  16.  Cluater  eCartlng  radius  vs.  position  at  vhlcb  it 
nucleates  (lower  curve)  and  corresponding  final  radius 
(upper  curve).  By  method  of  characteristics  with  Pg  •  20  torr 

0.1  cw  nossle  diameter. 
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Figure  18.  Clueter  initial  and  final  radiua  va. 
nucleation  position  for  Run  6  (evaporation  dieallowed) 
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Figure  20.  Cluster  initial  and  final  radiua  ya. 
nucleatlon  poaitlon  for  Run  7  (evaporation  diaallowed) 
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Figure  22.  Clueter  Initial  and  final  radiua  va. 
nucleatlon  poaition  for  Run  9  (evaporation  dlaallowed) 
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Another  drawback  of  our  approach  la  the  use  of  the  conventional  steady 
state  nucleatlon  theory.  Since  the  jet  expansion  Is  very  rapid •  It  Is 
probably  preferable  to  use  a  non-steady  theory,  following  the  evolution  of  the 
entire  spectrum  of  drop  sizes  and  solving  the  kinetic  equations  for 
condensation  and  evaporation  (Tang  &  Lu,  1986). 

He  have  also  considered  the  possibility  of  cluster  coalescence  to  explain 
the  presence  of  large  clusters.  This  explanation  appears  unlikely  since  the 
estimated  mean  free  path  for  clusters  at  Pq  *  20  torr  and  *  2140*’K  Is 
over  100  cm.  This  figure  was  arrived  at  by  considering  the  cross-sectional 
area  and  number  density  of  nuclei  when  the  mass  fraction  y  ~  0.01  as  In  case 
5  of  Table  I.  Even  If  y  ~  .1,  the  mean  free  path  would  be  much  too  long  to 
give  a  substantial  amount  of  coagulation  growth. 

Radiative  cooling  of  small  clusters  might  also  be  Important  In  determining 
their  temperatures,  and  therefore  their  ability  to  pass  off  excess  heat  of 
condensation.  To  our  knowledge,  this  possibility  has  not  been  considered 
elsewhere,  and  did  not  occur  to  us  until  the  modeling  work  was  largely 
completed.  However,  It  Is  certainly  plausible  that  the  high  temperatures  of 
the  particles  combined  with  their  high  surface  to  volume  ratio  could  lead  to 
some  significant  degree  of  cooling.  The  cooling  rate  due  to  thermal  radiation 
Is 


il  -  _  gAoT^ 
dt  me 

P 

where  e  is  the  emmlsslvlty,  A  the  surface  area,  a  the  Steffan-Boltxmann 
constant,  m  the  mass,  and  Cp  the  specific  heat.  This  can  be  Integrated  to 
give  the  temperature  as  a  function  of  time: 

T  -  (1/T  ^  +  9eOt/prc 
o  p 

As  an  example,  we  calculated  the  temperature  of  a  lOO  atom  aluminum 
cluster  (r  ■  0.8  nm)  which  Is  Initially  at  2100  K.  If  it  Is  assumed  that  e 
“  0.1,  then  the  cluster  will  cool  to  1798  K  In  10”^  seconds.  In  the  nozzle 
flow  of  an  expanding  jet  this  represents  about  2.5  mm  of  travel,  comparable  to 
the  size  of  the  condensation  zone.  Such  high  cooling  rates  could  not  only 
release  the  heat  of  condensation  but  may  act  as  a  sink  for  the  vapor  as  a 
whole,  increasing  the  superaturatlon  level. 

A  difficulty  arises  In  trying  to  apply  this  model  In  a  realistic  way 
because  the  emlsslvltles  are  unknown.  The  emmlsslvlty  of  a  metal  Is  closely 
related  to  Its  reflectance,  which  Is  In  turn  a  function  of  the  dielectric 
constant.  Ordinarily  metals  have  low  emmlsslvltles  (less  than  0.1)  and  are 
poor  radiators.  However,  Krelblg  and  Fragsteln  (1969)  modeled  the  optical 
properties  of  silver  particles  and  found  that  the  refractive  index  should 
Increase  dramatically  for  particles  smaller  than  the  mean  free  path  of 
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conduction  electrons  in  the  metal.  Using  their  model  we  calculated  that  a 
silver  sphere  of  1  nm  radius  will  have  a  reflectance  of  about  0.59,  as 
compared  with  0.96  for  bulk  silver,  at  a  wavelength  of  400  nm.  This  suggests 
an  emmlsslvlty  of  about  0.4,  leading  to  even  more  rapid  cooling  than  the 
example  of  the  paragraph  above.  Other  complicating  factors  such  as  wavelength 
and  temperature  dependence  and  relaxation  times  ought  to  be  considered  In  a 
complete  treatment  of  radiative  cooling,  but  our  estimates  show  that  it  Is  a 
potentially  Important  mechanism. 

4.5  Other  Mechanisms  for  Cluster  Growth 


An  alternative  mechanism  for  cluster  formation  and  growth  has  been 
proposed  by  Knauer  (1986,  1987).  It  postulates  that  nuclei  grow 
heterogeneously  on  Interior  surfaces  of  the  crucible,  then  detach  and  are 
entrained  In  the  vapor  flow  through  the  nozzle.  This  mechanism  has  two 
advantages  over  the  homogeneous  model:  first,  the  substrate  can  take  up  some 
or  all  of  the  heat  of  condensation  of  the  metal  vapor;  and  second.  It  should 
be  easier  for  atoms  to  nucleate  In  a  two-dimensional  space  than  In  three 
dimensions. 

In  Rnauer's  model,  surface  nucleatlon  and  growth  rates  are  calculated  for 
metal  vapor  which  arrives  at  a  high  flux  on  a  surface  slightly  cooler  than  the 
rest  of  the  crucible.  This  could  be  the  walls  of  the  nozzle  or  nearby 
interior  surfaces.  If  the  surface  Is  much  cooler  (more  than  15°C)  or  If  the 
clusters  form  strong  bonds  to  the  surface,  they  will  continue  to  grow  until  a 
continuous  layer  is  formed.  However,  If  the  bonding  Is  weak,  high  energy 
phonons  can  eject  clusters  In  balance  with  the  nucleatlon  rate.  A  number  of 
assumptions  regarding  the  energies  of  attachment  and  diffusion  of  atoms  and 
the  bonding  energy  to  the  surface  had  to  be  made  In  Rnauer’s  paper;  however, 
he  was  able  to  predict  that  silver  clusters  of  about  600  atoms  ought  to  be 
producible  In  reasonable  numbers  within  a  narrow  range  of  conditions. 

The  conditions  are  In  fact  very  narrow  and  difficult  to  establish  or 
measure.  Figure  23  shows  some  results  of  Rnauer's  calculations  which  Indicate 
how  closely  the  temperature  differences  must  be  controlled.  A  key  condition 
Is  that  the  metal  vapor  should  not  wet  the  crucible  material;  on  a  macroscopic 
scale  this  Is  easy  to  establish,  but  not  so  straightforward  for  aggregates  of 
a  few  tens  or  hundreds  of  atoms  on  a  rough  surface.  It  Is  also  quite 
difficult  to  measure  crucible  temperature  profiles  accurately  enough  In  the 
temperature  range  required. 

As  will  be  discussed  later,  we  performed  numerous  experiments  to  try  to 
establish  the  prescribed  conditions  and  did  not  find  evidence  of  clusters  In 
the  600  atom  range.  Another  group  working  at  Sohlo  (Windlschmann,  1987) 
obtained  a  similar  result.  However,  we  did  find  larger  particles  Issuing  from 
the  orifice,  and  were  able  to  correlate  their  presence  with  crucible  surface 
processes  by  studying  wall  deposits  on  rapidly  cooled  crucibles.  We  therefore 
regard  the  heterogeneous  growth  explanation  to  be  unproven  but  quite  plausible. 
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SURFACE  TEMPERATURE.  ^ 


Figure  23.  Cluster  size  vs.  crucible  teaperature 
differential,  froa  heterogeneous  theory 
(From  Knauer,  1987). 


5.0  EXPERIMENTAL  DEVELOPMENT 
5.1  Initial  Assumptions 

At  the  outset  of  the  Phase  I  investigation  we  acted  upon  several 
assumptions  which  guided  our  approach.  The  most  important  was  that  cluster 
formation  could  occur  in  a  pure  metal  vapor  expansion,  using  the  type  of 
nozzled  crucible  described  by  Takagi  (1975,  1984).  This  assumption  appeared 
reasonable  based  on  some  if  not  all  of  the  evidences  which  Takagi 's  group 
published  in  a  long  series  of  papers,  as  well  as  private  and  group  discussions 
held  with  members  of  his  research  team.  Credible  evidence  of  large  clusters 
has  recently  appeared  from  other  sources  (Knauer,  1987)  so  chat  this 
assumption  is  not  seriously  doubted. 
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k  second  assuaptlon  was  that  the  homogeneous  nucleatlon  and  growth 
Mchanism  proposed  by  Takagl  was  essentially  correct.  This  mechanism  is  known 
to  account  for  the  growth  of  water  and  rare  gas  clusters  in  supersonic 
expansions.  Also,  calculations  had  recently  been  published  (Usui,  1985)  which 
predicted  sufficiently  high  nucleation  rates  for  metal  vapors.  If  this  second 
assumption  were  true  then  it  should  have  been  possible  to  produce  clusters  in 
the  sise  range  which  has  been  reported;  Takagi's  papers  give  the  dimensions  of 
the  nozzles  and  the  crucible  stagnation  pressures  at  which  clusters  appeared. 
Surface  properties  would  then  play  no  part.  We  planned  to  make  our  major 
contribution  by  substantially  Improving  beam  characterization  techniques  and 
then  using  these  in  the  course  of  a  scientifically  based  thin  film  development 
program . 

Our  equipment  was  built  in-house,  tihich  allowed  us  to  quickly  bypass  many 
of  the  problems  which  other  U.S.  investigators  were  experiencing  with 
commercial  systems.  As  a  result  we  were  able  to  successfully  test  the  basic 
Takagl  crucible  configuration  using  a  well-characterized  analysis  and 
detection  system.  Because  of  the  very  low  upper  limit  which  was  then  placed 
on  cluster  population,  the  assumption  of  homogeneous  growth  was  set  aside  in 
favor  of  the  surface  growth  mechanism.  However,  we  were  then  faced  with  a 
greater  challenge  because  of  the  much  greater  difficulties  Involved  in 
Independently  measuring  and  controlling  a  multi-temperature  crucible.  The 
crucible  material  and  surface  condition  also  became  a  concern  at  that  point. 

A  third  assumption  was  that  clusters  in  the  supersonic  beam  should  have 
nearly  the  same  velocity  as  the  monomers.  This  implies  that  their  kinetic 
energies  will  be  proportional  to  their  masses,  so  that  energy  analysis  can  be 
used  to  obtain  a  mass  distribution.  At  the  temperatures  of  interest,  monomers 
accelerated  in  a  nozzle  expansion  will  have  an  average  translational  energy 
E  ■  5/2  kT  "  0.35eV.  If  N  “  1000  atoms  are  collectively  accelerated  before 
condensing,  the  cluster  should  have  E  *  5/2  HkT  *  350  eV.  This  method  is  not 
as  accurate  as  quadrupole  analysis  but  has  been  shown  to  give  good  agreement 
in  a  wide  range  of  experiments  with  rare  gas  clusters  and  seeded  beams.  It 
was  used  by  Yamada  for  silver  clusters;  typical  data  are  shown  in  Figure  2A. 

If  Instead  clusters  are  formed  on  crucible  surfaces  and  detach,  they  must  be 
accelerated  by  drag  forces  and  may  not  achieve  the  monomer  velocity.  However, 
we  have  performed  drag  calculations  which  show  that  1000  atom  clusters  should 
reach  97Z  of  the  monomer  velocity  while  traversing  a  1  mm  path  through  a 
nozzle. 

Arrangements  have  been  made  with  the  Microelectronics  Center  at  Iowa  State 
University  to  do  time-of-fll^t  mass  spectrometry  when  appropriate,  to  verify 
and  refine  the  mass  distributions  from  our  cluster  sources.  The 
tlme-of -flight  apparatus  was  not  completed  within  the  contract  period,  but  may 
be  useful  in  our  future  work. 


m-407/5-87 


37 


Figure  24.  Electrostatic  energy  analysis  of  silver 
cluster  beaa  (froa  Tawida,  1981). 


5.2  Initial  Hardware  Design 

A  ground*^p  approach  was  taken  in  developing  the  experimental  apparatus, 
Which  avoided  any  artifacts  or  difficulties  associated  with  ICB  equipment 
supplied  cooBBercially.  The  principal  investigator  also  had  the  benefit  of 
experience  in  design  of  ICB  components  from  earlier  work.  A  schematic  diagram 
of  our  source  apparatus  is  shown  in  Figure  25,  and  a  photograph  In  Figure  26. 

The  work  was  carried  out  in  an  Enerjet  UCV  18/6  diffusion  pumped  system, 
shown  in  Figure  27.  Special  care  was  taken  to  keep  chamber  pressures  below 
lOT^  torr  during  the  experiments  to  reduce  spurious  ion  signals  due  to 
background  gases.  A  residual  gas  analyser  was  used  to  monitor  actual  chamber 
gas  pressures  and  composition. 

Resistance  heating  was  used  for  the  crucibles.  This  was  much  simpler  than 
the  electron  bombardment  heating  used  by  Takagl,  and  kept  stray  electron  or 
Ion  currents  to  a  alniauu  in  the  chamber.  Initially,  boron  nitride  crucibles 
were  used  instead  of  graphite  because  BN  is  not  attacked  by  aluminum. 

Crucibles  were  machined  to  fit  in  commercially  available  crucible  beaters;  the 
lids  were  made  separately  and  In  early  versions  simply  rested  on  the  crucible 
body.  Extra  heat  shielding  was  added  to  Improve  temperature  uniformity.  A 
typical  crucible  Is  shown  In  Figure  28.  Thermocouple  and  pyrometer 
measurements,  and  modeling  of  radiative  heat  losses  from  the  lid,  determined 
that  the  lid  was  about  60^C  cooler  than  the  body  of  the  crucible.  The  lid 
orifice  had  a  1  mm  diameter  and  both  cylindrical  and  conical  profiles  were 
tried. 
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Figure  25. 


Schematic  diagram  of  apparatus. 
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Figure  27. 


High  vacuum  chamber  for  Phase  I 
cluster  beam  study. 
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Figure  28.  Boron  nitride  crucible  in  heater 


The  electron  Impact  ionizer  is  shown  partially  disassembled  in  Figure  29. 
It  Includes  the  grid  and  filament  from  a  conventional  ionization  gauge.  This 
made  construction  easy  and  had  the  advantage  that  its  ionization  efficiency  is 
known  or  can  be  calculated  for  various  gases  and  metal  vapors.  Because  the 
electrons  can  make  several  passes  through  the  beam  ionization  region  before 
being  intercepted  by  an  atom  or  collected  by  the  grid,  it  is  most  important 
that  ionizers  of  this  type  be  calibrated.  One  of  the  difficulties  in 
Interpreting  Takagl’s  results  is  that  no  such  calibration  is  ever  given,  only 
measurements  of  the  total  grid  current.  He  also  used  a  non-standard  grid 
configuration,  which  makes  it  virtually  impossible  to  even  guess  at  the  actual 
electron  current  densities.  The  characterization  of  our  ionizer  will  be 
discussed  in  Section  5.3.2. 
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Figure  29.  Ionization  stage  with  screen  removed 


A  key  requirement  which  we  placed  on  the  ionization  stage  was  that  it 
should  efficiently  transmit  all  species  in  the  beam.  Transmission  is  not  a 
problem  for  massive  clusters  which  obtain  roughly  300  eV  of  kinetic  energy 
from  the  nozzle  expansion.  However,  monomers  which  enter  the  ionizer  with  0.3 
eV  before  ionization  will  be  at  the  mercy  of  small  field  gradients  once  they 
are  charged.  Only  by  extracting  all  ions  could  we  hope  to  obtain  a 
quantitative  sampling  of  the  beam  constituents. 

To  meet  the  efficiency  requirement,  the  top  plate  of  the  ionizer  was 
electrically  Isolated  from  the  filament  and  grid  so  that  an  ion  extraction 
potential  could  be  placed  on  it.  Typically,  the  grid  was  maintained  at  75  to 
100  volts  and  the  ions  accelerated  toward  a  grounded  plate.  Extraction  of  an 
ionized  beam  was  problematic  at  first  with  the  open  cylindrical  grid.  The  ion 
currents  passed  to  the  Faraday  cup  were  too  low  and  the  energy  spread  too 
wide.  This  was  due  to  a  large  dead  spot  in  the  potential  gradient  within  the 
grid,  which  permitted  a  space  charge  to  build  up  and  inhibit  ion  flow.  The 
solution  was  to  weld  a  fine  nickel  mesh  around  the  outside  and  bottom  of  the 
grid  and  to  provide  a  small  grounded  tube  above,  to  extract  and  accelerate  the 
ions.  The  actual  geometry  was  modeled  by  solving  the  Laplace  equation  for  the 
potentials  within  the  grid  and  tube,  which  predicted  that  we  should  obtain 
mo  corner  ion  energies  of  65  to  75  eV.  This  was  confirmed  experimentally  as 
shown  in  Figure  30  for  ionized  silver. 
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Figure  30.  Energy  spread  of  extracted  ions 


Detection  of  the  beam  flux  and  current  was  done  with  a  crystal  rate 
monitor,  an  electrostatic  energy  analyzer,  and  a  shielded  Faraday  cup  (see 
Figure  25).  A  slot  was  cut  In  the  analyzer  and  the  contour  replaced  with  a 
fine  mesh  so  that  the  rate  monitor  could  measure  neutral  beam  flux  on-axls, 
while  the  analyzer  and  cup  measured  Ion  current  on-axls.  The  basic  principle 
of  a  cylindrical  electrostatic  analyzer  can  be  easily  understood  from 
Figure  31.  In  essence.  It  Is  a  device  which  creates  a  nearly  uniform  radial 
field;  by  setting  the  electrode  voltages  one  can  select  particles  of  a  given 
energy  while  others  collide  with  the  walls.  The  analyzer  was  constructed  In 
accordance  with  the  design  rules  by  Septler  (1967)  to  give  near-theoretical 
performance  and  focusing  along  the  silt  dimension.  The  analyzer  had  a  range 
of  0  to  800  eV  and  a  resolution  of  about  4Z. 

The  Faraday  cup  was  made  In  the  shape  of  an  elongated  box  to  match  the 
dimensions  of  the  analyzer  exit  slit.  This  maximized  the  detection  efficiency 
without  adding  much  to  Its  capacitance.  Shielding  was  provided  (not  shown  In 
Figure  26)  to  keep  out  stray  electron  currents.  Shielded  cables  routed  the 
current  signals  to  an  external  amplifier. 
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#  5.3.1  Cross  Section  for  Electron  Impact  Ionization 

An  Important  goal  of  this  work  was  to  obtain  quantitative  Information 
regarding  the  makeup  of  a  cluster  beam.  Electrostatic  methods  can  be  applied 
to  separate  charged  species  by  mass,  but  It  Is  much  more  difficult  to  analyze 
the  neutral  species.  One  method  which  has  sometimes  been  tried  Is  to  look  at 

#  TEM  grids  which  have  been  given  a  low-coverage  deposit  and  assume  that  the 
Islands  which  are  seen  correspond  to  actual  clusters.  This  method  has  worked 
well  for  gas-evaporated  particles  which  impact  at  low  velocity,  and  where  a 
continuum  of  sizes  from  microns  down  to  a  few  nanometers  can  be  produced  and 
studied.  We  were  hesitant  to  use  It  in  the  present  investigation,  however, 
because  we  expected  that  most  of  the  beam  flux  would  be  In  the  form  of  single 

#  atoms  rather  than  clusters.  Because  atoms  deposited  on  a  surface  will 
nucleate  Into  small  islands  of  their  own  it  would  be  difficult  to  draw  any 
firm  conclusions  from  such  studies. 


Instead,  one  can  try  to  find  the  neutral  population  by  counting  the 
cluster  ions  and  then  working  backwards.  This  requires  that  some  estimate  be 
#  obtained  for  the  ionization  cross  section  of  the  clusters,  and  also  that  the 

ionizing  apparatus  be  well  characterized.  This  Is  the  approach  which  we 
followed. 
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Few  If  any  data  are  available  for  the  ionisation  cross  sections  of 
clusters,  k  recent  and  coaprehensive  book  on  impact  ionisation  by  Mark  (1985) 
contained  only  two  or  three  pages  on  clusters «  and  made  no  mention  of  data  for 
metal  clusters.  In  the  absence  of  direct  data»  estimates  must  be  made  by 
geometric  considerations  and  extrapolations  from  atomic  ionisation  cross 
sections.  The  cross  section  for  ionisation  of  silver  atoms  is  shown  in 
Figure  32.  As  is  typical  for  virtually  all  atoms  and  molecules,  the  curve 
peaks  at  an  energy  several  times  the  ionisation  potential  and  then  gradually 
decreases  at  higher  energies.  Takagl  (1975)  published  a  curve  for  ion  current 
vs.  electron  energy  with  a  lead  cluster  beam  which  showed  the  same  general 
behavior;  other  than  his  curve  no  data  have  been  found  which  pertain  to  metal 
clusters . 

One  cross  section  estimate  can  be  had  from  the  geometrical  cross  section 
alone .  Thus 

o(N)  -  A  -  Tr(3Mv/4TT)^^^ 

where  A  is  the  frontal  area,  N  the  number  of  atoms  in  a  cluster,  and  v  the 
volume  per  atom.  Using  silver  clusters  of  1000  atoms  as  an  example  leads  to  a 
cross  section  of  8.1  x  10“^^  cm^.  This  method  was  used  by  Kulper  et  al. 

(1981)  in  their  study  of  cluster  sources.  It  is  obviously  a  crude 
approximation,  since  the  cross  section  can  have  no  energy  dependence. 

Another  estimate  (and  probably  an  upper  bound)  comes  from  treating  the 
atoms  in  a  cluster  as  completely  independent  particles.  Then  a(N)  ■ 

Na(l).  Again  using  lOOO  atom  silver  clusters,  this  leads  to  a  ■  4.6  x 
l(rl3  cm2  at  75  eV. 


Figure  32.  Electron  impact  ionisation  cross 
section  for  silver  (from  Crawford,  1967). 
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A  third  way  which  seens  preferable  is  to  treat  the  atoms  as  independent, 
but  to  include  a  penetration  depth  factor.  Mean  free  path  data  for  electrons 
in  solids  can  be  used  such  as  those  given  in  Somarjai  (1981).  For  75  eV 
electrons  in  silver  this  depth  is  0.5  nm,  considerably  less  than  the  expected 
cluster  diameters  of  about  3  nm.  Multiplying  the  depth  by  the  surface  area  of 
the  cluster  and  the  atom  density  gives  the  number  of  atoms  sampled  by  the 
electron,  about  230  for  a  1000  atom  cluster.  Then  a  ■  230  x  a(l)  ■  1.1  x 
10"^^  cvr .  Interestingly  enough,  this  value  is  quite  close  to  that  given 
by  the  frontal  area.  However,  this  holds  true  only  near  the  energy  at  which 
the  atomic  ionization  peaks.  At  higher  or  lower  energies  this  method  will 
reflect  the  general  tendencies  of  the  atomic  ionization  curve,  as  modified  by 
the  penetration  depth  factor.  At  350  eV  this  gives  o(l000)  ■  5.6  x  10“^^ 
cm^  for  the  cluster. 

A  valid  objection  which  might  be  raised  is  that  the  metal  clusters  really 
ought  to  behave  more  like  a  bulk  metal.  Ue  considered  several  ways  of 
accounting  for  this,  principally  by  applying  classical  and  empirical  models 
for  atomic  ionization.  The  atomic  models  were  modified  by  replacing  the 
atomic  ionization  potential  with  the  metal  work  function  and  then  summing  over 
the  quasi-independent  atoms.  As  a  test,  we  tried  to  reproduce  the  energy 
dependence  which  was  shown  in  by  Takagl  for  lead  clusters.  With  all  the 
models  tried,  the  peak  occurred  at  much  too  low  an  energy.  Whether  this  is 
due  to  the  weaknesses  of  the  model  or  the  available  data  is  uncertain. 

We  expect  that  better  experimental  data  and  models  will  someday  be 
available  for  cluster  cross  sections.  At  present,  there  does  not  seem  to  be 
any  better  way  of  estimating  these  quantities  than  by  the  method  of  summing 
over  Independent  atoms  within  the  penetration  depth  of  the  electrons. 

5.3.2  Detection  Efficiency  of  Apparatus 

Careful  measurements  of  the  performance  of  the  ionizer  and  electrostatic 
analyzer  were  needed  so  that  we  could  make  quantitative  statements  about  the 
beam  composition.  These  measurements  included:  1)  total  rate  of  gas  ion 
production  vs.  chamber  pressure;  2)  total  rate  of  metal  ion  production  vs. 
flux  or  deposition  rate;  and  3)  the  fraction  of  ions  which  were  extracted  from 
the  ionizer,  passed  the  entrance  slit  of  the  analyzer,  and  were  transmitted  to 
the  Faraday  cup. 

The  rate  of  gas  ion  production  was  important  both  as  a  check  on  operation 
and  so  that  such  currents  could  be  subtracted  from  the  metal  ion  current  later 
on.  This  was  done  by  temporarily  placing  a  biased  collector  wire  in  the 
center,  just  as  in  a  normal  ion  gauge.  After  pumping  the  chamber  down  to  less 
than  10~'  torr  and  outgassing  the  ionizer,  the  chamber  valve  was  closed  and 
the  ion  collector  current  read  as  the  pressure  slowly  increased.  Linear 
relationships  were  established  between  collector  current  and  pressure,  and 
between  collector  current  and  grid  current,  with  the  result  that  1-23 
mA/torr-mA.  The  residual  gas  was  primarily  water  vapor,  which  is  nearly 
equivalent  to  nitrogen  in  ionization  probability.  For  comparison,  the 
published  ionization  gauge  sensitivity  is  10  mA/torr-mA  for  nitrogen.  This  is 
reasonable  agreement  considering  that  element  spacings  can  be  quite  critical 
in  gauge  tubes. 
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Total  aetal  Ion  production  was  measured  in  a  similar  way,  except  that  beam 
flux  replaced  residual  pressure  as  a  parameter.  The  result  was  that  1  ■  13.2 
nA/mA-I/s  for  silver,  where  the  flux  was  measured  by  a  crystal  rate  monitor  17 
cm  above  the  source.  In  a  typical  beam  experiment  with  a  5  £/s  rate  and  10  mk 
of  grid  current  the  metal  Ions  would  account  for  660  nA.  This  is  well  above 
the  roughly  20  nA  contribution  from  Ionized  residual  gas  at  10~'  torr. 

With  the  collector  wire  removed,  the  Ionizer  extraction  and  analyzer 

transmission  efficiencies  were  measured.  For  the  metal  Ions  and  a  0.5  by  2  cm 

• 

slit,  42Z  of  the  beam  or  5.6  nA/mA-A/s  was  collected.  This  is  much  higher 
than  the  IZ  or  less  typical  of  most  mass  spectrometer  Ionizers,  In  which  the 
extraction  field  gradient  must  be  reduced  to  maintain  high  resolution. 

The  evidence  which  was  sought  for  cluster  formation  was  an  Ion  current 
signal  at  a  substantially  higher  energy  than  the  75  eV  or  so  given  to  all 
charged  species  In  the  Ion  extraction  process.  A  cluster  %rlth  N  ■  1000  will 
have  a  kinetic  energy  of  about  1000  x  5/2  kT  ■  350  eV  before  being  Ionized, 
and  425  eV  after  Ionization  and  extraction.  In  order  to  find  the  mass 
fraction  F(N)  of  the  beam  which  Is  In  the  form  of  N-atom  clusters,  the  beam 
current  I(N)  at  high  energy  can  be  compared  with  the  low  energy  monomer 
current  1(1)  provided  that  the  Ionization  cross  sections  are  taken  Into 
account : 


F(W)  ,  0(1)  KN)  „ 

fHT  iTTT  *  " 

where  a(N)  Is  the  cluster  cross  section,  ct(1)  the  monomer  cross  section, 
and  the  factor  of  N  accounts  for  the  mass  per  Ion.  It  Is  assumed  that  the 
monomer  fraction  F(l)  accounts  for  most  of  the  beam  mass  (F(l)  =  1). 


The  detection  limit  for  cluster  population  depends  on  the  cluster  size  and 
also  on  the  signal  to  noise  ratio  at  the  Faraday  cup  amplifier.  With  the 
analyzer  tuned  away  from  the  monomer  ion  peak  the  noise  limit  was  0.006  nA 
peak-to-peak ,  five  orders  of  magnitude  below  the  typical  metal  Ion  current. 

The  finite  resolution  of  the  energy  analyzer  must  also  be  taken  Into  account 
by  calculating  the  equivalent  mass  range  sampled  within  the  FWHM.  For  N  ■ 
1000,  dM  “  40.  Thus  the  beam  fraction  within  this  range  will  be  given  by 


F( 980-1020) 


0(1)  1(980-1020) 

o(1000)  *  1(1) 


X  1000 


Inserting  o(l000)  -  1.1  x  10”^^,  o(l)  -  4.6  x  10“^^,  1(980-1020)  -  6 
pA  (the  noise  limit)  and  1(1)  ■  450  nA  (typical)  gives  F(980-1020)  ■  5.6  x 
10^^.  A  similar  calculation  for  N  ■  500,  ct(500)  ■  6.9  x  10“^^,  gives 
F(488-512)  ■  4.4  X  10~^.  To  a  good  approximation,  these  results  can  be 
extrapolated  to  find  the  detectable  mass  fractions  over  broader  size  ranges: 

Fdet^ ^00-600)  =  3.7  x  10“^ 

Fdet(®00-1200)  =  5.6  X  10“^ 
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The  800  to  1200  atoB  range  was  chosen  because  It  represents  a  fair  sample  of 
the  cluster  distribution  reported  by  Tsaada  (1981).  The  400  to  600  range 
corresponds  to  the  results  of  Knauer  (1987).  These  two  examples  and  the  text 
above  show  what  we  believe  our  detection  limits  were  during  the  noeele  beam 
experiments  and  how  we  arrived  at  them. 

5.4  Performance  of  Noes led  Crucibles 


5.4.1  Results  of  Cluster  Detection  Experiments 

Data  published  by  Takagl  show  that  clusters  are  most  numerous  at  a  vapor 
pressure  of  a  few  torr,  and  virtually  disappear  outside  the  range  1  to  5 
torr.  Vapor  pressure  tables  for  liquid  metals  show  considerable  scatter,  but 
a  3  torr  pressure  ought  to  occur  in  the  range  1820  to  1900  K  for  aluminum  and 
1575  to  1675  K  for  silver.  The  vapor  pressure  changes  by  about  a  factor  of 
three  with  a  100  K  change  for  either  material.  In  most  of  our  experiments  two 
tungsten/ rhenium  thermocouples  were  used,  one  in  the  body  of  the  crucible  and 
the  other  in  the  noczled  lid. 

A  more  direct  means  of  determining  the  vapor  pressure  is  through  measuring 
the  rate  m  of  mass  flow  from  the  crucible  and  using  the  following  relation  for 
a  sonic  nozzle: 


m  -  ✓  aY(Y<-l)/2RTg  Pg  (2/( 4^ 

where  A^  is  the  nozzle  cross  sectional  area,  a  the  molar  mass,  R  the  gas 
constant,  and  P^  and  Tq  the  pressure  and  temperature  in  the  crucible,  y 
is  taken  to  be  5/3,  the  specific  heat  ratio  for  a  monatomic  gas. 
Experimentally,  the  mass  flow  rate  was  determined  to  an  accuracy  of  lOX  or 
better  by  weighing  the  metal  charge  in  the  crucible,  evaporating  part  of  it  at 
a  nearly  constant  rate,  and  then  ire-weighlng  afterwards.  As  an  example,  a 
flow  rate  for  silver  of  1.6  x  10~^  gm/s  was  obtained  with  a  1.25  mm  orifice, 
which  by  the  above  relation  gives  Pq  ■  4.8  torr.  A  thermocouple  on  the 
crucible  body  read  Tq  ■  1706  +  20  K,  well  within  the  range  expected  from 
vapor  pressure  tables. 

In  our  first  series  of  experiments,  aluminum  was  used  because  of  its 
usefulness  as  an  optical  and  electronic  material.  It  was  evaporated  from 
boron  nitride  crucibles  to  avoid  reactions,  which  other  cluster  beam  workers 
had  noted  when  using  aluminum  in  graphite.  Electrostatic  analyzer  data  were 
taken  from  0  to  800  eV  for  PqCAI)  *  0.5,  1.6,  and  3.5  torr.  A  10  mA  grid 
current  at  75  eV  was  used  for  ionization,  giving  a  monomer  peak  of  over  a 
hundred  nA  at  3.5  torr.  No  other  peaks  were  seen.  The  lid  in  this  case  was 
about  50^  cooler  than  the  body  of  the  crucible,  which  should  not  be 
Important  if  clustering  occurs  in  the  gas  flow  rather  than  on  surfaces. 
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In  the  absence  of  a  detectable  peak*  an  upper  Halt  for  cluster  population 
can  be  established.  This  Halt  Is  deteralned  by:  1)  the  signal  to  noise  ratio 
of  the  detection  electronics;  and  2)  the  Ionisation  cross  section  of  a 
cluster.  Both  of  these  were  discussed  In  the  preceding  sections  and  a 
detection  Halt  arrived  at  for  silver.  In  the  early  alualnua  ezperlaents  the 
noise  level  was  soaewhat  higher,  about  0.1  nA,  but  a  slallar  analysis  leads  to 
the  following  upper  Halts  for  cluster  aass  fractions  In  the  500  and  1000  atoa 
ranges: 

Pai( 400-600)  <  6.1  X  10"^ 

Pai(800-1200)  <9.3  x  10“^ 

After  these  Initial  results  we  received  a  preprint  outlining  the 
heterogeneous  growth  aechanlsa  (Knauer,  1987).  This  aodel  was  discussed  In 
Section  4.5.  Aaong  Its  predictions  Is  that  alualnua  should  not  fora  clusters 
In  a  cluster  beaa  source  because  alualnua  will  wet  virtually  any  crucible 
aaterlal  at  high  teaperatures .  This  laplles  that  clusters  growing 
heterogeneously  on  crucible  surfaces  will  adhere  too  strongly  to  be  ejected  by 
phonon  excitation.  This  led  us  to  try  silver  Instead  of  alualnua  In 
subsequent  experlaents . 

In  the  silver  experlaents,  special  care  was  taken  to  keep  the  crucible  Hd 
and  orifice  only  slightly  cooler  than  the  crucible  body  teaperature.  This  was 
done  to  keep  conditions  within  the  range  predicted  by  Knauer  for  heterogeneous 
nucleatlon.  This  degree  of  teaperature  control  Is  not  easy  to  achieve  at  high 
teaperatures  because  of  the  large  radiation  losses  and  the  difficulty  of 
aalntalnlng  good  themal  conduction  paths,  plus  the  difficulty  of  asking  good 
teaperature  deteralnatlons.  If  the  heterogeneous  ;r*^ichanlsa  Is  correct,  the 
best  approach  to  crucible  design  would  Include  a  separate  heater  and  power 
supply  for  the  lid  and  orifice  region.  This  was  not  possible  In  the  Phase  I 
study.  However,  by  using  two  heat  shields  over  the  crucible  lid  and  a  long 
threaded  connection  between  the  parts  the  lid  teaperature  could  be  kept  to  -20 
+  lO^C  of  the  body.  A  finite  difference  aodel  of  the  radiant  and  conductive 
heat  transfer  was  used  to  predict  the  Hd  teaperature;  this  served  as  an  aid 
to  design  and  a  check  on  the  reasonableness  of  the  theraocouples.  An  exaaple 
of  the  crucible  designs  Is  shown  In  Figure  33. 

Crucibles  of  both  boron  nitride  and  a  carbon- loaded  graphite  were  used  In 
the  silver  beam  experlaents.  Also,  some  of  the  noszles  had  a  simple  1  x  1  am 
cylindrical  profile  while  others  had  a  tapered  section  as  was  shown  In 
Figure  33.  In  the  course  of  checking  various  Ideas  many  trials  were  made, 
which  eventually  covered  most  of  the  pressure  range  between  0.3  and  16  torr. 
Figures  34  and  35  are  relatively  unadorned  original  data  typical  of  the  runs 
made,  showing  control  over  temperature  and  deposition  rate,  and  an 
energy-analyzed  beam  current  spectrum. 
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Figure  33.  Crucible  design  £or  minimizing 
temperature  differential. 


Considerable  evidence  was  accumulated  which  pointed  toward  the 
heterogeneous  mechanism  for  cluster  growth,  including  the  ejection  of  fine 
droplets  under  certain  conditions.  These  results  will  be  discussed  In  the 
next  section.  However,  the  electrostatic  analyzer  data  did  not  show  the 
presence  of  clusters  In  the  nanometer  size  range  of  about  20  to  2000  atoms 
each. 

An  upper  bound  for  the  silver  cluster  population  la  then  given  by  the 
detection  limit,  which  was  derived  In  the  previous  section.  Thus 

P^(  400-600)  <3.7  X  10"^ 

P.  (800-1200)  <  5.6  X  10"^ 

Ag  — 

when  pQ  •  9  torr.  At  Pq  ■  3  torr  the  upper  limit  is  about  a  factor  of  two 
higher  because  of  the  lower  monomer  current. 
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Figure  34.  Crucible  teaperatures  and  deposition 
rates  during  typical  run. 
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Faraday  Cup  Current 


Figure  35. 


Example  of  electrostatic  energy 
analyser  results. 
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The  flow  field  from  the  nossle  had  been  modeled  previously,  as  described 
In  Section  3.2.  From  this  Information  It  Is  possible  to  predict  the  thickness 
distribution  of  a  film  deposited  on  a  substrate.  As  a  test  of  the 
reasonableness  of  the  gas  flow  modeling,  several  semitransparent  silver  films 
were  deposited  for  comparison.  Film  thickness  profiles  were  determined  from 
densitometry  and  calculations  using  the  optical  constants  of  silver.  The 
results  for  nossle  beam  deposition  with  "  1  torr  are  shown  In  Table  II. 

The  agreement  Is  found  to  be  quite  good,  with  the  predicted  profile  being 
peaked  In  comparison  with  an  effusive  source  and  the  measured  profile  peaked 
slightly  more.  The  calculations  were  carried  out  up  to  an  angle  where  more 
than  half  the  beam  mass  was  accounted  for  (M0  >  0.5),  although  the 
substrate  subtended  only  20*. 


Table  II. 


Predicted  Versus  Measured  Thickness  Profile 


6f inoi 

Actual 

Predicted 

Effusive 

R/RN 

fd^T 

Me 

(nm) 

(nm) 

Source  (nm) 

0.0 

0.00 

0.00 

50.0 

50.0 

50.0 

0.1 

4.04 

0.01 

49.1 

49.3 

49.5 

0.2 

7.98 

0.04 

46.7 

47.6 

48.1 

0.3 

11.96 

0.09 

42.1 

44.4 

45.8 

0.4 

16.11 

0.16 

38.9 

40.2 

42.6 

0.5 

20.61 

0.25 

34.7 

38.4 

0.6 

25.34 

0.36 

28.3 

33.4 

0.75 

33.40 

0.54 

17.6 

24.3 

r/Rk  - 
®flnal  ■ 


Me  - 


Starting  position  of  streamtube  In  nossle,  expressed  as  radial 
distance  from  centerline.  (R/R]i]  ~  0.0  on  centerline,  *  1.0  at 
nossle  wall). 

Final  angle  of  streamtube  from  R,  calculated  by  method  of 
characteristics . 

Fraction  of  total  beam  mass  within  ^flnal. 
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5.4.2  Evidence  for  Heterogeneous  Growth  of  Fine  Particles 

The  experiments  described  In  the  previous  section  seem  to  show  that 
significant  metal  cluster  growth  cannot  occur  through  homogeneous  nucleatlon 
mechanisms.  In  spite  of  this  result,  other  mechanisms  cannot  be  ruled  out. 

One  which  seems  plausible  Is  heterogeneous  nucleatlon  at  Inner  surfaces  of  the 
crucible.  Knauer  (1987)  has  worked  out  a  reasonable  theory  for  this  mechanism 
which  can,  given  certain  assumptions,  predict  that  clusters  form  (see  Section 
4.5).  The  requirements  Include  a  slightly  cooler  surface  which  experiences  a 
high  flux  of  metal  vapor  and  non-wetting  (weak  bonding)  of  cluster  material  to 
the  crucible.  If  the  temperature  difference  Is  small,  clusters  should  detach 
from  the  walls  before  growing  appreciably  and  may  then  be  entrained  In  the 
vapor  flow. 

The  theory  gives  temperature  ranges  In  which  clusters  could  be  formed  by 
this  mechanism.  Several  of  the  cluster  detection  experiments  were  run  In  or 
near  this  range,  with  /ST  *  -20  +  10*C.  While  no  clusters  of  N  less  than 
2000  atoms  were  observed.  It  should  be  added  that  the  predicted  range  Is 
probably  somewhat  sensitive  to  the  data  and  assumptions  fed  Into  the  model, 
and  only  the  general  trends  of  the  theory  have  been  verified  by  Knauer  so  far. 
Also,  the  ultimate  cluster/particle  size  should  be  a  sensitive  function  of  the 
crucible  temperature  differential,  as  was  shown  In  Figure  23.  Besides 
electrostatic  analyzer  measurements,  some  TEM  grids  were  briefly  exposed  to 
the  beam  and  analyzed  by  Ren  Westmacott  at  the  National  Electron  Microscopy 
Center  at  La%rrence  Berkeley  Laboratory.  These  did  not  show  evidence  of 
clusters. 

Positive  evidence  for  the  surface  (heterogeneous)  growth  mechanism  was 
obtained  by  observations  of  larger  particles.  In  some  runs  droplets  were 
ejected  from  the  orifice,  ranging  In  size  from  lOO  nm  to  several  microns.  A 
silver  coating  which  Includes  many  such  particles  Is  shown  In  Figure  36.  Some 
of  the  smallest  particles  seen  are  shown  In  Figure  37,  in  which  a  film  of  3  nm 
average  thickness  was  deposited  on  carbon.  The  many  small  patches  In  the 
background  are  Islands  of  silver  which  nucleated  on  the  substrate  from  the 
monomer  component  of  the  beam.  It  Is  clear  that  the  particles  are  not  merely 
the  upper  end  of  the  Island  size  distribution. 

The  ejected  particles  could  be  directly  correlated  with  the  metal  deposits 
formed  on  the  underside  of  the  crucible  lid  near  the  orifice.  The  low  thermal 
mass  of  our  crucibles  allowed  us  to  cool  them  at  nearly  l5°C/s,  rapidly 
lowering  the  metal  evaporation  rate  and  preserving  most  of  the  deposits. 
Surface  properties  appeared  to  correlate  with  particle  presence  In  two  ways : 

1)  material  combination,  and  2)  surface  finish  of  the  orifice  region. 

The  material  combinations  Investigated  were;  aluminum  In  BN  crucibles; 
silver  In  BN;  and  silver  In  carbon-filled  graphite.  The  aluminum  both  wet  and 
reacted  with  BN,  as  shown  in  Figure  38,  and  no  particles  were  ejected  from  the 
orifice.  Likewise,  no  particles  were  seen  with  silver  In  BN,  even  though  this 
Is  considered  to  be  a  non-wetting  combination.  Figure  39  of  the  orifice 
region  shows  that  particles  do  frrm  there  but  that  the  area  of  surface  contact 
Is  still  quite  large.  However,  a  much  higher  contact  angle  of  about  135*  was 
observed  for  silver  on  carbon  (Figure  40)  and  particles  were  ejected  under 
certain  conditions.  This  correlation  of  particle  ejection  with  low  surface 
bonding  strongly  suggests  that  Knauer 's  theory  of  heterogeneous  cluster  growth 
Is  basically  correct. 
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Figure  36. 


Silver  coating  formed  from  nozzle  beam  source. 
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Figure  37.  Silver  particles  and  Island  structure  from  nozzle  beam  source. 
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Figure  38.  Alunloun  reacted  with  boron  nitride  crucible  lid. 
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Figure  41.  Effect  of  surface  finish  on 
droplet  formation  in  crucible. 


Another  Interesting  observation  was  that  particles  were  much  more  likely 
to  be  ejected  when  the  orifice  region  had  a  relatively  coarse  finish. 

Figure  41  shows  the  underside  of  a  crucible  lid  in  which  the  center  portion 
has  been  smoothed.  A  few  large  particles  appear  in  the  smoothed  portion  but 
it  is  otherwise  clear.  This  is  hard  to  explain,  since  roughuess  ought  to 
#  bring  about  more  complete  wetting  and  stronger  bonding.  A  reasonable 

explanation  is  that  roughness  holds  the  droplets  more  tenaciously,  during  the 
growth  period,  so  that  they  reach  a  size  which  is  observable  with  our  SEM. 

The  very  large  particles  are  probably  the  result  of  coalescence.  In  any  case, 
the  correlation  of  particle  ejection  with  surface  properties  supports  the 
heterogeneous  growth  hypothesis. 

5.4.3  Growth  of  Ultrafine  Particles  or  Clusters  by  Gas  Quenching 

The  apparent  difficulty  of  making  metal  clusters  with  a  Takagl-t3rpe  source 
led  us  to  explore  other  methods.  This  was  a  natural  extension  of  the  original 
scope  of  work  and  has  already  opened  up  new  options  for  future  research. 
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Gas  quenching  Is  one  technique  for  obtaining  particles  reliably  and  with 
considerable  control  over  the  final  particle  size.  In  Its  simplest  form>  a 
metal  can  be  heated  In  an  open  boat  to  a  vapor  pressure  of  a  few  torr  and  the 
vapor  quenched  In  a  quasl-statlc  Inert  gas.  The  nucleatlon  and  growth  of 
particles  Is  then  controlled  by  the  pressures  and  by  the  velocity  at  which  the 
gas  passes  over  the  meral  source.  Typically »  convection  currents  In  the  meter 
per  second  range  are  allowed  to  control  the  gas  velocity  and  the  particles 
drift  toward  collection  surfaces.  A  source  of  this  type  Is  shown  In 
Figure  42.  Granqvlst  (1976)  performed  a  thorough  study  of  the  pressure 
dependences  and  found  that  more  or  less  monodlsperse  particles  could  be  formed 
with  any  desired  size  from  a  few  nanometers  to  hundreds  of  nanometers. 

This  simple  form  Is  similar  In  principal  to  the  pulsed  laser  vaporization 
sources  now  used  widely  In  cluster  science.  The  pulsed  sources  provide  tiny 
quantities  of  small  clusters  (N  £  40) ;  the  clusters  and  carrier  gas  can  then 
be  formed  Into  supersonic  beams  for  spectroscopic  studies.  Velocities  of 
hundreds  of  meters  per  second  are  obtained. 

An  ultraflne  particle  source  suitable  for  coating  deposition  might  be 
developed  by  using  parts  of  both  methods.  A  boat  or  oven  can  efficiently 
vaporize  large  quantities  of  material  Into  a  quenching  gas,  and  the  resulting 
particles  entrained  In  a  nozzle  flow  to  provide  Initial  acceleration. 

Figure  43  Illustrates  this  concept.  A  number  of  Issues  arise  In  properly 
designing  such  a  source,  such  as:  maintaining  the  desired  gas  velocity  over 
the  metal  source  (to  minimize  or  control  convective  forces);  keeping  the 
particles  In  the  center  of  a  laminar  flow  to  avoid  orifice  fouling;  creating 
flow  convergence  patterns  Into  the  nozzle  which  do  not  lead  to  centrifugal 
settling  of  particles  out  of  the  flow;  differential  pumping  of  the  downstream 
portion  to  remove  or  recycle  the  carrier  gas  and  allow  deposition  under 
moderate  vacuum.  The  proposal  for  Phase  II  of  this  program  addresses  these  In 
some  detail.  Ultimately,  a  combination  of  gas  flow  modeling  and  experimental 
development  will  be  needed. 

To  conclude  the  Phase  I  development  effort,  a  proof  test  was  conducted  of 
gas-quenched  particle  growth.  Silver  was  evaporated  from  our  nozzle  beam 
source  at  a  high  rate  Into  1  torr  of  argon.  An  open  crucible  could  have  been 
used  to  the  same  end  (see  Figure  42),  but  the  nozzled  source  gave  the  metal 
vapor  a  more  uniform  flow  velocity  In  the  nucleatlon  zone  and  should  therefore 
result  In  a  narrower  distribution  of  particle  sizes.  Of  course,  beyond  the 
first  few  millimeters  the  partlcle/gas  mixture  should  have  slowed  to  a  speed 
determined  by  convection  In  the  static  chamber. 

An  SOI  photo  of  the  resulting  particles  Is  shown  In  Figure  44.  The 
deposit  consisted  of  fine  particles  of  approximately  50  to  100  nm,  loosely 
agglomerated  as  expected.  Formation  of  a  dense  coating  will  require 
entrainment  of  the  particles  In  a  converging  nozzle  flow  along  the  lines  of 
Figure  43.  However,  the  basic  Idea  of  gas-quenched  growth  seems  to  be  sound, 
and  the  nozzled  crucible  for  the  metal  vapor  may  be  useful  In  developing  the 
overall  source. 
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For  optical  aaterlal  applications  a  saaller  particle  aay  be  aore 
desirable.  The  literature  indicates  we  could  reach  this  range  with  a  aetal 
source  which  gives  a  lower  super saturation,  such  as  an  ICB-type  crucible  with 
a  larger  orifice  at  lower  teaperature.  A  lower  argon  pressure  or  hellua 
instead  of  argon  should  also  reduce  the  particle  else  (Granqvlst,  1976). 

Gas-quenched  techniques  are  especially  interesting  because  of  the  range  of 
particle  slses  obtainable.  Thus,  it  aay  be  possible  to  enoraously  broaden  the 
range  and  possibilities  of  "cluster"  deposition.  At  the  2  -  3  na  end  of  the 
scale  are  large  clusters  of  roughly  a  thousand  atoas,  slallar  to  those 
considered  by  Takagi  for  epitaxial  deposition.  With  a  high  fraction  of 
surface  atoas,  it  is  intuitively  plausible  that  such  clusters  should  break  up 
into  individual  atoas  on  iapactlng  a  substrate  at  high  velocity  (see  Figure 
2).  Larger  clusters  of  5  -  10  na  or  aore  alght  aore  properly  be  teraed 
ulrafine  particles.  With  fewer  surface  atoas,  these  should  retain  aore  of 
their  particulate  character  and  aay  fora  coatings  by  flattening  or  "splattlng" 
on  a  substrate  (see  Figure  3).  This  would  be  wore  in  line  with  the  structure 
of  plasaa  sprayed  coatings,  but  on  a  auch  finer  scale.  Because  plasaa  sprayed 
coatings  do  not  show  the  columnar  structure  of  aost  thin  fllas,  deposition  of 
ultraflne  particles  aay  provide  significant  advantages  for  optical  filas. 
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Figure  44.  SEM  photo  of  50-100  nn  silver  particles 
produced  with  nozzled  crucible  and  argon  quenching* 
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Personnel 

Principal  Investigator 

Arther  C.  Day,  M.S.  Applied  Physics,  University  of  Washington,  1983. 
Thesis:  "Analysis  of  Surfaces  by  Rutherford  Backscattering.” 

Contributors 


Julian  A.  Domaradzkl ,  Ph.D.,  Physics,  Institute  of  Theoretical 
Physics,  Warsaw,  Poland. 

6.  S.  Knoke,  Ph.D.,  Aeronautics  and  Astronautics,  University  of 
Washington,  1975.  Thesis:  "Comhusted  Ludwieg  Tubes  for  Gaadynamlc 
Laser  Research." 

R.  P.  Turcotte,  Ph.D.,  Physical  Chemistry,  Arizona  State  University, 
1969.  Thesis:  "1.  High-Temperature  X-Ray  Studies  of  Widely 
Nonstolchlometrlc  Regions  in  the  Praseodymium  Ozlde-Ozygen  System; 

11.  Structure  and  Thermal  Stability  of  the  Rare  Earth  Carbonates  and 
Formates." 

Inventions 


Section  4.4.3  of  this  report  discussed  two  concepts  that  may  be  considered 
inventions.  Their  development  was  outside  the  original  scope  of  work  and 
did  not  result  in  additional  contract  change;  however,  they  are  a  logical 
extension  of  the  Phase  1  effort  and  are  relevant  to  future  work  as 
contemplated  In  the  Phase  II  proposal.  They  are: 

1.  Gas-quench  assisted  foirmatlon  of  ultraflne  prarticles  with  nozzle 
beam  source. 

2.  Laminar  flow  enclosure  for  transport  of  ultraflne  particles  and 
separation  from  quench  gas. 

Publications 


The  following  papers  cover  work  performed  on  or  related  to  the  Phase  1 
contract  and  are  planned  for  submission  In  August  1987 : 

"Emlsslvities  and  Cooling  Rates  for  Hot  Clusters,"  A.  C.  Day,  to  be 
submitted  to  Applied  Physics  Letters  or  Journal  of  Applied  Physics. 

"Mechanisms  for  Cluster  Formation  and  Deposition:  Limitations  and 
Possibilities,"  A.  C.  Day,  J.  A*  Domaradzkl,  G.  S.  Knoke,  to  be 
submitted  to  Journal  of  Vacuum  Science  and  Technology. 
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